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While current nuclear parton distribution functions (nPDFs) from global fits to experimental 
data are spatially homogeneous, many experimental observables in nucleus-nucleus collisions are 
presented in terms of centrality cuts. These cuts can be related to impact parameter using the 
Glauber theory and it is thus usual in the description of such observables to convolute an assumed 
impact parameter distribution with the homogeneous nPDFs. In this study we use the Gribov 
theory of nuclear shadowing supplemented with information from diffraction to model the impact 
parameter distributions of nuclear shadowing ratio in the small- a; region. The modeled distributions 
are applied to the description of the centrality dependence of observables in deuteron-gold (d-|-Au) 
collisions at y^sjviv = 200 AGeV. 
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I. INTRODUCTION 

Nuclear parton distribution functions (nPDFs) encode 
the modifications to free nucleon parton distributions 
due to the complex, many-body environment in the nu- 
cleus. They are essential in the application of perturba- 
tive Quantum Chromodynamics (pQCD) to the descrip- 
tion of many relevant observables in relativistic nucleus- 
nucleus collisions. Most of the currently available nPDFs 
have been obtained from global fits to diverse exper- 
imental data. These nPDFs from global fits are homoge- 
neous (independent of impact parameter), and are thus 
functions of the Bjorken variable x and squared momen- 
tum transferred Q^. An exception is [6] which is based on 
the Gribov theory Q and available in both homogeneous 
and inhomogeneous forms. 

Experimental data on nucleus-nucleus collisions are of- 
ten presented in terms of centrality classes. These classes 
can be mapped into corresponding intervals in impact pa- 
rameter space using the Glauber formalism. In order to 
adequately describe these data, inhomogeneous (impact- 
parameter-dependent) parton distributions are required. 
Due to the inherent constraints in the determination of 
nPDFs, especially in the gluon sector, it is pragmatic 
not to be limited in the choice of nPDF sets. Thus in 
order to take advantage of the different sets available, 
there is need for phenomenologically robust and realis- 
tic prescriptions for generating inhomogeneous nPDFs. 
This is usually carried out by convolution of impact pa- 
rameter distributions with homogeneous nPDFs. One 
such prescription is given in [8] where two variants of 
impact parameter distributions are considered. The first 
variant assumes the distribution is proportional to the 
local nuclear density (generally taken as a Wood-Saxon 
nuclear density). The second takes the distribution as 
proportional to the nuclear path length (Glauber thick- 



ness function). A simpler distribution is used in Q where 
a hard sphere dependence is assumed. 

Our objective in this study is to model impact param- 
eter distributions which can be used with the homoge- 
neous nPDFs in the spirit of d The general idea is 
to determine the impact parameter distributions which 
give a good description of the nuclear shadowing ratio in 
the region of small Bjorken x. These distributions are 
then applied in conjunction with a homogeneous nPDF 
set to calculate the centrality dependence of some observ- 
ables in deuteron-gold collisions. A major simplification 
underlying the application is the assumption that the dis- 
tributions derived from small- a: shadowing are applicable 
to modifications at all values of x. As in [ll|, [l^] we 
employ a generalized form of the Gribov theory, incor- 
porating the real part of the diffractive scattering ampli- 
tude. Experimental data for nuclear shadowing at small 
x are available from the NMC[ll,[il| and E665gl[il| ex- 
periments. These experimental data are all at small Q^: 
thus the use of information from diffractive scattering of 
real photons (Q^ = o GeV^ at FNAL ^) and quasi-real 
photons (g2 < 0.01 GeV^ at HERA [18,]) afl?ords a good 
approximation. 

In order to present a reasonably self-contained analy- 
sis, we include all the relevant details from our previous 
treatments. The paper is thus organized as follows: in 
Sec. In] we review the basic formalism of Gribov theory 
as applied to shadowing in the small Bjorken- a; regime. 
This section is essentially the same as in ■ We present 
the results of our calculations for the shadowing ratio and 
the impact parameter distributions in Sec. lIIII In Sec. lIVI 
we illustrate a typical application of our results by using 
the distributions to describe the centrality dependence 
of the nuclear modification factor, pseudorapidity asym- 
metry, and hadron ratio in ultrarelativistic deuteron-gold 
(d-|-Au) collisions. We conclude in Sec. |Vl 
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II. NUCLEAR SHADOWING FROM 
GENERALIZED GRIBOV THEORY 

A. Nuclear Shadowing Ratio 

In nuclei, for small values of the Bjorken variable x 
{x < 0.1), the nuclear structure functions are sig- 
nificantly reduced compared to the product of the mass 
number A and the free nucleon structure function F2 ■ 
Since the virtual photon-nucleus cross section is propor- 
tional to F2 , then equivalently the virtual photon- nucleus 
cross section is less than A times the one for free nucleons, 
a^-A < Aa^*M- This phenomenon is generally known as 
nuclear shadowing in the strict sense. Similar behavior 
is observed for real photons at sufficiently high energies 
(i' > 3 GeV). Thus the nuclear shadowing ratio, defined 
as F2 /{A * F2) or alternatively as a^*A/{A * ct^'at), is 
less than unity. 

The (virtual) photon-nucleus cross section is separable 
into a part which accounts for the incoherent scatter- 
ing from individual nucleons, and a correction (shadow- 
ing correction) from the coherent interaction with several 
nucleons: 



'■7* A 



Z cF^tp + {A — Z) fT^.„ + Saj*A 



(1) 



The single scattering part is the incoherent sum of 
photon-nucleon cross sections, where Z is the nuclear 
charge number, and cr^-p and a-y*n are the photon-proton 
and photon-neutron cross sections, respectively. The 
multiple scattering correction is expressible as an expan- 
sion in the number of nucleons in the target involved in 
the coherent scattering (n > 2). The dominant contri- 
bution to nuclear shadowing comes from double scatter- 
ing, since the probability that the propagating hadronic 
excitation coherently interacts with several nucleons de- 
creases with the number of nucleons. 



We define the shadowing ratio as 



{A - Z) a. 



7*n 



(5(7, 



Z (J^'p + [A — Z) a 



(2) 



7^^ n 



We neglect the small difference between the photon- 
proton cross section a^p and the photon-neutron cross 
section (T^„ in this study, and thus denote by a^- n the 
generic photon-nucleon cross section. The shadowing ra- 
tio can thus be written as 



n^A = 1 



(5(7. 



7* A 



Act. 



7* AT 



(3) 



Thus, the evaluation of the shadowing correction, Sa-y^A, 
is central to the calculation of the shadowing ratio. In the 
next section we utilize the Gribov theory in a generalized 
form to determine Sa^-^A- 




FIG. 1: Double scattering contribution to Sa^ 



B. Shadowing Correction From Generalized 
Gribov Theory 

The original formulation of Gribov is generalized by 
including the real part of the diffractive scattering am- 
plitude. We denote by 77 the ratio of the real to imagi- 
nary parts of the diffractive scattering amplitude. In this 
generalized form the shadowing correction at the level of 
double scattering is given by 



(5(7, 



7* A 



A{A-1) 



2^2 



167^7^e 



dzi / dz2 



p^^\b, zi; 6, Z2) exp i i 



dM 



l-f 772 



'7*Ar 



^ dM^dt 

jzi - Z2) 

X 



(4) 



with a^i]^ the photon-nucleon diffractive cross section. 
The coherence length. A, is given by A = 2z//M^ for 
real photons. As illustrated in Fig. [1] a diffractive state 
with invariant mass Mx is produced in the interaction 
of the photon with a nucleon located at position (b, zi) 
in the target. The hadronic excitation is assumed to 
propagate at fixed impact parameter b and to interact 
with a second nucleon at Z2. The probability to find two 
nucleons in the target at the same impact parameter is 
described by the two-body density p\ {b, zi;b, Z2) nor- 
malized as / d'^r d'^r' p'"J\f,f') = A^. The phase factor, 
exp{i[(zi — Z2)/X]} in Eq. ^ implies that only difi^rac- 
tively excited hadrons with a longitudinal propagation 
length larger than the average nucleon-nucleon distance 
in the target, A > li ~ 2fm, can contribute significantly 
to double scattering. The limits of integration define the 
kinematically permitted range of diffractive excitations, 
with their invariant mass Mx above the two-pion produc- 
tion threshold and limited by the center-of-mass energy 
W = y/s of the scattering process. 

We approximate the two-body density Pa {b, zi; b, Z2) 
by a product of one-body densities p^(r)/9^(r') since 
short-range nucleon-nucleon correlations are relevant in 
nuclei only when Z2 — 21 is comparable to the range of the 
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short-range repulsive part of the nucleon-nucleon force, 
i.e. for distances ^ 0.4 fm. However, shadowing is neg- 
hgible in this case and therefore short-range correlations 
are not important in the shadowing domain. 

With increasing photon energies or decreasing x down 
to X <C 0.1, the longitudinal propagation length of diffr ac- 
tively excited hadrons rises and eventually reaches nu- 
clear dimensions. Thus for heavy nuclei, interactions of 
the excited hadronic state with several nucleons in the 
target become important and should be accounted for. 
Following we introduce an attenuation factor with an 
effective hadron-nucleon cross section, (loff- The shadow- 
ing correction can thus be written as 



Sa^'A = — — ^ IGnTZe 



2^2 



77^ 



dzi 



dz2 



4,„2 ^ dMjrdt 



p^Hb, zi;b, Z2) exp 



A 



exp 



-{1/2){1 ~ iri)acs dzpA{b,z) 



(5) 



The effective hadron-nucleon cross section, (Toff in Eq. (O 
is defined as 



167r 



fcff 



dM^ ^ ^^'"^ 



^ dM^dt 



(6) 



where g^n is the photon-nucleon cross section. The de- 
tails of this approach and the approximations inherent in 
the definition of cFcS are treated thoroughly in Q. For 
vector mesons as the intermediate hadronic excitations, 
we take ay n as cTcff in the attenuation factor in Eq. ([5]) , 
where avN is the vector meson-nucleon scattering cross 
section. 



C. Impact Parameter Profile 

The shadowing correction ([5]) can be expressed as 



6a^,A = JdHr^ib) 



(7) 



where 



87^7^e 



(1 - z??)2 
1 -I- 772 

j2 diff 



dz-i / dzo / dMy — 



py{b,zi;b,Z2) exp 



jzi - Z2) 

X 



exp |-(1/2)(1 - i77)crcfi y dzpA{b,z] 



(8) 



We make use of the fact that the shadowing correction, 
Sa^tA^ is a real number and write 



27rbJ-|(6) 

(5tT. 



JdbGiib)^ 1 



(9) 



which is reminiscent of the thickness function distribu- 
tion: 



J^^ 2TTbJdzpA{b,z) ^ _^ 



(10) 



The integrands in Eq. (|9]) and Eq. (fTO|) for both light and 
heavy nuclei are shown in Sec. IIIII 



D. Diffractive Dissociation 

The shadowing correction, (Jcr^.^i, which is the major 
component in the evaluation of the shadowing ratio and 
impact parameter profile, involves the differential diffrac- 
tive dissociation cross section. The determination of this 
cross section is the subject of this section. 

The analysis by the HI Collaboration [l^ divides the 
HERA diffractive dissociation data into effectively three 
intervals in M^: the first interval (0.16 - 1.58) GeV2 
contains the region of the low-mass vector mesons (p, lu 
and </)), the second interval (1.58 — 4.0 GeV2) covers the 
p' resonance region, while the third interval (M|- > 4.0 
GeV2) is that of the high- mass continuum. Due to the 
similarity in the theoretical description of the first and 
second intervals, we combine them together. The differ- 
ential diffractive cross section is thus written as a sum 
over contributions from two mass intervals. 



dMj.dt 



E 



tP=iO V—p,uj,<p,p' 



dMj^dt 



fjfjcont 



t«0 



dMjrdt 



(11) 



In the following, we briefly summarize the various ap- 
proximations applied in these intervals. 



1. Low-mass vector mesons and the p' resonances 

We utilize the generalized vector meson dominance 
(VMD) [lH to describe the contribution of the low-mass 
vector mesons to the differential diffractive cross section: 



dM\dt 



IGtt 



Ml 



x) 2 



(12) 



with H^(M|.) the vector meson part of the photon spec- 
tral function H(AfD, which is given by 



h(a4) 



1 ^7(e^ 



hadrons) 



127r2 (7(e+e — > p'^ p ) 



(13) 



In Eq. avN is the vector meson-nucleon cross sec- 

tion and e2/47r = 1/137 is the fine structure constant. 
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The w and 4> mesons are narrow and thus well approxi- 
mated by delta functions. Their contribution to the pho- 
ton spectral function can be written as 



X) 



mv 
9v 



5{M 



where my and gv, {V — are the mass and the 

coupling constant of the uj and </) mesons, respectively. 

The p-meson, unlike the uj and </> mesons, has a large 
width due to its strong couplin g to two-pion states. We 
have followed the approach in [20| and taken this into 
account through the tt+tt" part of the photon spectral 



function: 



) = j^e(Mj-4„e)(i-M) 



3/2 



If, (m 



X) 



(15) 



where ttItt is the mass of the pion and Mx — M^^-^ is 
the invariant mass of the tt+tt" pair. The pion form 
factor, is taken from Ref. [2l|. A full discussion is 
given in [2^. We compared the result from the delta 
function approximation to this more exact calculation 
and found that taking into account the width of the p- 
meson increases the differential diffractive cross section 
by about 10%. 

The p' resonance region contains the p(1450) and 
p(1700) mesons which were formerly classified as the 
p(1600) [13]. The FNAL data show an enhancement in 
this region and we treat this enhancement in terms of an 
average p' resonance, corresponding to the earlier clas- 
sification of p(1600), as done in Ref. [l3|- We use the 
available information on the p(1600) from Ref. [l^ in a 
VMD-type calculation to evaluate the contribution from 
this region. The average p' resonance should have a finite 
width, but encouraged by the fact that a delta function 
in the case of the p gives a good approximation to the 
full- width result, we employ a narrow resonance approx- 
imation for the p(1600): 



dMjrdt 



tsiO 



16n Mj. 



'VN 



with 



nV(Mi) = (^) 5{Ml-mi.) 
9v 



(16) 



(17) 



and V = p(1600). 



2. High-mass continuum 

A full treatment of both the FNAL data and the 
HERA data in this region has been carried out by 
the HI Collaboration in Ref. [l^, using the triple- 
Regge model of photon dissociation. Here we consider 



only the aspects relevant to the present study. There 
are two diffractive terms in the triple-Regge expansion: 
the triple-pomeron (PPP) and the pomeron-pomeron- 
reggeon (PPIR) terms. The subleading reggeons have 
the quantum numbers of the p, w, 02 and /2 mesons and 
their trajectories are approximately degenerate. They 
are generally referred to as p, w, a and / mesons and their 

isospin, signature and C- and G- parities are p(l h), 

a;(0 ), a(l + H — ) and /(0 + ++). The pomeron has 

P(0 -I- -I— 1-) and is thus identical to the / meson, leading 
to interference. We neglect such interference effect in the 
present study. 

The differential dissociation cross section can thus be 
written as 



dV 
dM2 dt 



Gppp(O) 



Gip]p]R,(0) 



2Nj2ajp(0)-2 ^_B(W^M^)^ 



(18) 



where B{W'^,Ml) = 2hp^ + 2a'^ \n{W'^/Ml) and hpw is 
the proton-pomeron slope parameter. a'_^ is the slope 
of the pomeron trajectory while and are the 
pomeron and (effective) reggeon intercept respectively. 
We use the values in [l^ for these parameters. Note 
that the value of a^(0) (ajp(O) = 1.068 ± 0.0492) agrees 
within error with the soft pomeron intercept in Ref. (23j 
(q!^ (0) ~ 1.081). The usual triple-pomeron approxima- 
tion corresponds to putting Gppir(O) = in Eq. (fT8t . 
The contributions from the subleading reggeons are small 
at both NMC and E665 energies [12[, thus we use the 
triple-pomeron approximation in this study. We use the 
value of Gip]p]p(0) — 9.0 from the analysis in jl2 |. 



III. 



RESULT: SHADOWING RATIO AND 
IMPACT PARAMETER PROFILE 



The treatment outlined in the last section is now ap- 
plied to calculate the shadowing ratio and the impact 
parameter profile. The basic expression is Eq. (jSj, which 
involves the ratio of the real to imaginary amplitudes 77, 
the photon- nucleon cross section a-yN^ the nuclear den- 
sity Pa, and the effective cross section in terms of 
the diffractive dissociation cross section. 

We use the e nerg y-independent 77's for the vector 
mesons from Ref. 19]. For both p and w mesons, rj takes 
values between and —0.3. Here we take rjp = rj^ = —0.2 
in accordance with Ref. [l^. The results of our calcula- 
tion are not very sensitive to the precise values of ?7p((^) . 
For the meson, we take rj^ — 0.13 24]. Due to lack 
of information, we take r]p(isoo) = 0. For the high- mass 
continuum, we follow Ref. [6| and define rj^p as 



:K(0)-l) 



(19) 



using the result of Gribov and Migdal [25j . 

We use the Donnachie-Landshoff parameterization of 
(J |23| as the generic photon- nucleon cross section 
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For the nuclear densities three-parameter Fermi {3pF) 
distributions are apphed: 

P^"^^ ^P" l + e(r-RA)/d ' (20) 

with the parameter values taken from Ref . \2d\ . For mass 
numbers ^ < 20 a harmonic oscillator (HO) density dis- 
tribution may be more appropriate than the SpF dis- 
tribution. For uniformity, we use the 3pF distributions 
for the whole mass range in light of the fact that un- 
certainties associated with other parameters are at least 
comparable. 



B. Impact Parameter Profile At Small x 

We use Eq. ([SJ to determine the impact parameter 
profile for nuclei in the mass range 3 < A < 239 and at 
small X (a; ~ 10"''). In Fig. |3]we present the profile as 
a function of impact parameter b for four representative 
nuclei: '^c, 63Cu, '^^^u, and ^ospb. These results are 
shown in the figure as solid curves. 

As noted in the introduction, two prescriptions are 
given in |8J]: a local density dependence and a Glauber 
thickness function scaling. At small Bjorken-x the thick- 
ness function dependence is probably more physical. Its 
impact parameter distribution, from Eq. (|10p . is shown 
as dashed curves in Fig. |31 For both the shadowing cor- 



A. Shadowing Ratio At Small x 



At very small a; (a; ~ 10 ^) the E665 experiment has 
four data points: ^^c, 40Ca, ^^^Xe, and ^ospb. The re- 
sults of our calculation and the experimental data are 
shown in Fig. [21 For small A the shadowing ratio de- 
creases rapidly with A, while for large A the decrease is 
more gradual. The result is in good agreement with ex- 
periment, and reproduces quite accurately the trend of 
the small x shadowing with mass number A. 
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FIG. 2: (Color Online) Shadowing ratio at small x. The 
dashed curve is the result of our calculation. The shaded 
triangles are data from the E665 experiment. 
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FIG. 3: (Color Online) Impact parameter profile at small x. 
The solid and dashed curves are the results from the present 
study and thickness function dependence respectively. 

rection and thickness function, the profiles have narrow 
widths and sharp peaks for light nuclei, with the widths 
and peaks broadening out as the mass number increases. 
Since the profiles are normalized to unity, it is easy to see 
their differences: the profiles from the shadowing correc- 
tion have taller peaks and a steeper dropoff to zero than 
those from the thickness function dependence. 

Fig. [3] only shows the profiles as a function of b for four 
nuclei. In Fig. 2] we show the full result of our calculation 
of the profile from shadowing correction as a function 
of both the mass number A and impact parameter b. 
The features noted in the preceding paragraph are clearly 
manifested. 
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A. Collinear factorization 




FIG. 4: (Color On line) Impact parameter profile at small x 
as a function of mass number A and impact parameter b. 



The cross section for the d+Au — > h+X reaction, with 
respect to pseudorapidity rj and transverse momentum 
Pt can be written as 



dr](PpT 



J d^b (fs dxadxbdzc td{s) tAui\b- s\) 

abed 

Fa/d{Xa,Q'^,S,z) Fh/Au{xb,Q'^,\b- s|,z') 

da{ab->-cd) Dh/c{zc,Q}) 



dt 



-sS{s + t + u) (23) 



where Xa and Xb are parton momentum fractions in 
deuteron and gold, respectively, and Zc is the frac- 
tion of the parton momentum carried by the final-state 
hadron h. The factorization and fragmentation scales are 
Q and Qf respectively. Here 



(24) 



tAis) = dzpAis,z) 



is the Glauber thickness function of nucleus A, with the 
nuclear density distribution, pa{s, z) subject to the nor- 
malization condition 



d^s dzpA{s, z) = A 



(25) 



For applications to deuteron-gold collisions (Sec. ITV)) 
it is convenient to define the dimensionless form of the 
impact parameter distribution. A suitable form for the 
shadowing correction (sc) is 



-^I(O) 



(21) 



while for the thickness function (tf) we can define 



JdzpA{b,z) 
JdzpA{0,z) 



(22) 



These forms are not unique; any suitable value of b can 
be used in the denominator to render the profiles dimen- 
sionless. 



IV. APPLICATION TO DEUTERON-GOLD 
(D+AU) COLLISIONS 

We now apply the impact parameter profile discussed 
in Sec. IIIIBI to calculate the centrality dependence of 
some observable quantities in deuteron-gold collisions: 
nuclear modification factor for 7r°, R^au^ pseudorapidity 
asymmetry, Yas^: for pions and protons, and hadron-to- 
pion ratio. 



The quantity da[ab — > cd)/dt in Eq. (j23p represents 
the perturbatively calculable partonic cross section, and 
Dh/ci^cQ f'^) stands for the fragmentation function of 
parton c to produce hadron h, evaluated at momentum 
fraction Zc and fragmentation scale Q / . 

Since we are interested in the description of centrality- 
selected data using nPDFs which are spatially homoge- 
neous, we need to incorporate some form of impact pa- 
rameter dependency. We follow the approach in [1, 
and write 

Fa/A{x, Q'^,s, z) = fa/N{x, Q^)S{A, X, Q'^ ,s, z) . (26) 

Here fa/N{x, Q^) is the nucleon parton distribution func- 
tion, which can be expressed as 

fa/N{x.Q^) = |/a/p(a;,g2) + (l-|)/a/„(x,g2) , (27) 

with fa/p{x,Q'^) {fa/nix, Q"^)) being the proton (neu- 
tron) parton distribution function as a function of 
Bjorken x and factorization scale Q. The inhomogeneous 
shadowing function is defined as 

S{A,x,Q^,s,z) = 1 + N^[S'{A,x,Q^) - l]$(s,z) (28) 

with $(s, z) dimensionless and a suitable normaliza- 
tion constant such that 

J Jd^sdzpA(s, z)SiA, X, Q^,s, z) = S\A, x, Q^) (29) 
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Here S'{A,x,Q'^) denotes the homogeneous shadow- 
ing function. Note that the homogeneous nPDFs are 
the product of the homogeneous shadowing function 
S'{A,x,Q^) and the nucleonic PDFs fa/N{x,Q^). The 
dimensionless function $(s, z) furnishes the required im- 
pact parameter dependency. We use the profile derived 
from the shadowing correction, W^ib) (Eq. (HI])), sub- 
ject to the normalization condition in Eq. (j29p . 

The density distribution of the deuteron is ob- 
tained from the Hulthen wave function j27| (as in 
Ref. [15), while the density distribution for gold is 
taken as a Woods-Saxon distribution with parameters 
from Ref. [1^. We fix the scales as Q = Qf — 
Pt, where pr is the transverse momentum of the fi- 
nal hadron. The partonic differential cross sections, 
da{ab — cd)/dt were evaluated at leading order (LO). 
We note that, if a if factor was used to approximate 
the effects of higher orders, these efi'ects would cancel in 
the ratios calculated in the present study. We use the 
Eskola-Paukkunen-Salgado (EPS08)[2] shadowing func- 
tions with the MRST2001 leading order (LO) parton dis- 
tribution functions (PDFs) [291. For the fragmentation 
functions we use the DSS set [s^, IMj for pions and pro- 
tons. 

As stated earlier experimental data are often presented 
in terms of centrality classes, the most central class corre- 
sponding to head-on collisions. The Glauber description 
of heavy ion collisions afford a connection between cen- 
trality classes and impact parameter cuts. An analytic 
expression valid to a very good approximation is given 



^dAu 



(30) 



ratio of the d-|-Au and proton-proton (pp) cross sec- 
tions, normalized by the average number of binary colli- 
sions, (Nbin). 

The PHENIX Collaboration has presented results for 
the nuclear modification factor for neutral pions [s^ in 
the pseudorapidity interval |t;| < 0.35 and four centrality 
classes: - 20%, 20 - 40%, 40 - 60%, and 60 - 88%. We 
use these classes and the results of our calculations are 
displayed in Fig. [SJ together with the experimental data. 




10 40 



p,[GeV] 



where c is the centrality, b the impact parameter, and 
'^dAu is inelastic deuteron-gold cross section. 



B. Nuclear Modification Factor 

The d-|-Au nuclear modification factor for a hadron h, 
R^Au^ is defined as 

[Nbin) dr] d^pT I arj d^pr 

where the average number of binary collisions, (Nbin) in 
the various impact-parameter bins is given by 

{Nmu) = {<Jnn TdAuib)) . (32) 

Here cr™^ is the inelastic nucleon-nucleon cross section 
and 

TdAu{b) = jd^s td{s) tAuiSb- si) (33) 

represents the deuteron-gold nuclear overlap function. 
The nuclear modification factor -R^^„ is thus just the 



FIG. 5: (Color Online) Nuclear modification factor RdAu for 
neutral pion in different centrality classes. Filled triangles 
denote experimental data from the PHENIX Collaboration 

Our calculation overshoots the data for the central 
class (0 — 20%). This is in part a reflection of the 
fact that the minimum bias EPS08 overpredicts the data 
by about 5%. The data shows the influence of shad- 
owing (i?2^,j < 1) for Pt < 2.5 GeV, antishadowing 
> 1) for 2.5 < PT < 6 GeV, and the EMC ef- 
fect (i?5^„ < 1) for PT > 6 GeV. While our calculation 
follows this general trend, the data shows a narrower win- 
dow for the antishadowing peak relative to calculation. 
Away from the central class, the antishadowing peak pro- 
gressively broadens with increasing centrality, and this is 
adequately reproduced by calculation. The agreement 
with data is in fact quite reasonable for the other cen- 
trality classes under consideration. 
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C. Pseudorapidity Asymmetry 

As the mechanisms for hadron production in d+Au 
coUisions may be different at forward rapidities (deuteron 
side) and backward rapidities (gold side) , it is of interest 
to study ratios of particle yields between a given rapidity 
value and its negative in these collisions. The STAR 
Collaboration has measured pseudorapidity asymmetries, 
defined as 



YAsy = , (Au-side) 



(d-side) , (34) 



in d+Au collisions for several identified hadron species 
and total charged hadrons in the pseudorapidity inter- 
vals \-q\ < 0.5 and 0.5 < \r]\ < 1.0 Here we con- 
sider the pseudorapidity asymmetry for both pions and 
protons in different centrality classes for the two STAR 
pseudorapidity intervals. 

Fig.|6]shows the result of our calculation of the pseudo- 
rapidity asymmetry for pions (upper panels) and protons 
(lower panels) for the interval \ri\ < 0.5 and centralities 
— 20% and 40 — 100% respectively. For the pions, our 



data very well for the — 20% (central) class. It is off by 
a few percent for the peripheral class. 

Two general remarks are in order: first, the magnitude 
of the pseudorapidity asymmetry is rather small in this 
interval, generally less than or of the order of 10%. This 
is manifested both by experiment and calculation, and 
is due to the fact that the considered interval is close to 
midrapidity. Second, the calculation exhibits a greater 
degree of asymmetry for the central class relative to the 
peripheral. 

The result of our calculation for the pseudorapidity in- 
terval 0.5 < \ri\ < 1.0 is presented in Fig. [T] Away from 
midrapidity, the asymmetry is appreciably enhanced and 
this is reflected by both experiment and calculation. 
Also, as observed in the previous interval, the magnitude 
of the asymmetry decreases with increasing centrality. 
Overall, the agreement with data is manifestly good for 
both pions and protons and in both central and periph- 
eral classes. 



asy 0.8 




2 4 6 8 10 120 2 4 6 8 10 12 

p,[GeV] 



FIG. 6: (Color Online) Centrality dependence of pseudora- 
pidity asymmetry for pions (upper panels) and protons (lower 
panels) in the interval \ri\ < 0.5. Filled triangles denote ex- 
perimental data from the STAR Collaboration jsj). 

calculations are in good agreement with data for both 
centrality classes, especially for the 40 — 100% (periph- 
eral) class where the agreement is manifestly excellent. 
In the case of protons, the calculation reproduces the 



3=^0.8 




2 4 6 8 10 120 2 4 6 8 10 12 

p,[GeV] 



FIG. 7: (Color Online) Centrality dependence of pseudora- 
pidity asymmetry for pions (upper panels) and protons (lower 
panels) in the interval 0.5 < \ri\ < 1.0. Filled triangles denote 
experimental data from the STAR Collaboration [341 ]. 



D. Hadron Ratio 

As a final example of the application of the impact pa- 
rameter profile to d-|-Au collisions, we consider the cen- 
trality dependence of the ratio of hadrons to pions, as 
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measured by the PHENIX Collaboration [351. We de- 



note this ratio as "R-^^^^ and write 



(h + +h-)/2 



dAu 



(35) 



The centrality classes are the same as for the 7r° nuclear 
modification factor, viz: - 20%, 20 - 40%, 40 - 60%, 
and 60 — 88%. Fig. [5] shows the result of our calculation 
for these four centrality classes. The trend also mirrors 



that observed for R 



dAv 



Our calculations overshoot the 



R 



3 
2.5 
2 
1.5 

1 

h/jt 0-5 
dAu 

2.5 

2 
1.5 

1 

0.5 







M| ,11 

0-20% J 

+ffj — 1 — 


:J0^ 


20-40% ^ 

Hi\ \ — \-\-+ 


1 — 


T 40-60% : 


\ — I-H+ 


60-88% _ 




r s ^ 




t ~~ 

i ^ 

T s 




1,1 , , , 






1 


10 40 

Pt 


1 

[GeV] 


10 40 



i3 
-2.5 
■2 



0.5 



2.5 

2 



0.5 





FIG. 8: (Color Online) Centrality dependence of hadron-to- 
pion ratio. Filled squares denote experimental data from the 
PHENIX Collaboration fsl. 



3 < ^ < 240 and for small a; (a; ~ 10~^) using a suit- 
ably generalized form of the Gribov theory. The photon 
diffractive dissociation cross section, which is the main 
input to our calculation, has been parameterized as a 
function of the invariant mass of the difFractively pro- 
duced hadronic excitation in two effective mass intervals. 
For the low-mass vector mesons and p' resonances, vector 
meson dominance model is used while the high-mass con- 
tinuum is treated within the framework of triple-Regge 
theory. Relevant model parameters are taken from ex- 
periments and earlier studies. 

The calculated result for the shadowing ratio agrees 
nicely with small x shadowing data from the E665 exper- 
iment and also reproduces quite well the trend of shad- 
owing with mass number A. We compare the resultant 
impact parameter profiles for four representative nuclei 
with the Glauber thickness function dependence which 
has been previously applied in the literature. We also 
present the profiles as a function of both mass number A 
and impact parameter b for 3 < A < 240 and < 6 < 12 
fm. 

We illustrate the utility of the impact parameter dis- 
tributions by applying them to calculate the centrality 
dependence of three different observables in deuteron- 
gold collisions: the tt" nuclear modification factor 
pseudorapidity asymmetry YAsy for pious and protons, 
and the hadron-to-pion ratio, 'R-'^au- Except for the over- 
shoot at the most central class, our calculation gives an 
adequate description of the centrality behavior of both 

-^dAu ^'^d 'R-'Iau ^ measured by the PHENIX Collabora- 
tion. The agreement with the STAR data for the central- 
ity dependence of pseudorapidity asymmetry is also good 
for both pions and protons. This is encouraging consid- 
ering the fact that the distributions have been derived 
from the small- a; shadowing ratio, and thus in a strict 
sense reflects the behavior in the shadowing region, i.e 
X < 0.1. 



data for the most central class (0 — 20%) especially at 
higher px- The agreement with data is more reasonable 
for the other centrality classes. 



CONCLUSION 
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